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the hadron masses of the non-strange vector mesons, the hyperons and the nucleon 
in those nuclei. We find a new, simple scaling relation for the changes of the hadron 
masses, which can be described in terms of the number of non-strange quarks in the 
hadron and the value of the scalar mean- field in a nucleus. 
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I INTRODUCTION 

One of the most exciting topics in nuclear physics is the study of the variation of hadron 
properties as the nuclear environment changes. In particular, the medium modification 
of the light vector-mesons is receiving a lot of attention, both theoretically and experi- 
mentally. Recent experiments from the HELIOS-3 and the CERES |0| collaborations 
at the SPS/CERN energies have shown that there exists a large excess of the e^e~ pairs 
in central S + Au collisions. Those experimental results may give a hint of some change 
of hadron properties in nuclei [^. Forthcoming, ultra- relativistic heavy-ion experiments 
(eg. at RHIC) are also expected to give significant information on the strong interaction 
(QCD), through the detection of changes in hadronic properties (for a review, see Ref.[^). 

Theoretically, lattice QCD simulations may eventually give the most reliable infor- 
mation on the density and/or temperature dependence of hadron properties in matter. 
However, current simulations have been performed only for finite temperature systems 
with zero baryon density 0] . Therefore, many authors have studied the hadron masses in 
matter using effective theories: the vector dominance model QCD sum rules and 



the Walecka model P, [TTI , and have reported that the mass decreases in the nuclear 
medium (see also Ref. [|T^] ). 

In the approach based on QCD sum rules, the reduction of the mass is mainly due to 
the four-quark condensates and one of the twist-2 condensates. However, it has been sug- 
gested that there may be considerable, intrinsic uncertainty in the standard assumptions 
underlying the QCD sum-rule analyses [1^. In hadronic models, like Quantum hadrody- 
namics (QHD) the on-shell properties of the scalar (a) and vector {uo) meson with 
vacuum polarization were first studied by Saito, Maruyama and Soutome [§], and later by 
many authors P, |lO|, |rT|. (Good physical arguments concerning the uj meson in medium 
were found in Ref.[|l^.) The main reason for the reduction in masses in QHD is the 
polarization of the Dirac sea, where the anti-nucleons in matter play a crucial role. From 
the point of view of the quark model, however, the strong excitation of N-N pairs in 



medium is difficult to understand |IT5| . 
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Recently Guichon, Saito, Rodionov and Thomas |T6| have developed an entirely dif- 
ferent model for both nuclear matter and finite nuclei, in which quarks in non-overlapping 
nucleon bags interact self- consistently with (structureless) scalar [a) and vector [uj and 
p) mesons in the mean-field approximation (MFA) - the quark-meson coupling (QMC) 
model. (The original idea was proposed by Guichon in 1988 [0. Several interesting 



applications to the properties of nuclear matter and finite nuclei are also given in a se- 



ries of papers by Saito and Thomas ||T5|, |T9|, ^.) This model was recently used to 
calculate detailed properties of static, closed shell nuclei from ^^O to ^°^Pb, where it was 
shown that the model can reproduce fairly well the observed charge density distributions. 



neutron density distributions etc. |2T|]. Blunden and Miller ||22| have also considered a 
model for finite nuclei along this line. 

To investigate the properties of hadrons, particularly the changes in their masses in 
nuclear medium, one must also consider the structure of the mesons, as well as the nucleon. 



Saito and Thomas pOf have studied variations of hadron masses and matter properties 
in infinite nuclear matter, in which the vector-mesons are also described by bags, but 
the scalar- meson mass is kept constant, and have shown the decrease of the hadron mass. 
Now it would be most desirable to extend this picture to finite nuclei to study the changes 
of hadron properties in the medium, quantitatively. 

Our main aim in this paper is to give an effective Lagrangian density for finite nuclei, 
in which the structure effects of the mesons {a, uo and p) as well as the nucleon are 
involved, and to study quantitative changes in the hadron (including the hyperon) masses 
by solving relativistic Hartree equations for spherical nuclei derived from the Lagrangian 
density. (Using this model, we also calculate some static properties of closed-shell nuclei.) 
In the present model the change in the hadron mass can be described by a simple formula, 
which is expressed in terms of the number of non-strange quarks and the value of the scalar 
mean-field (see also Ref.||20|). This is accurate over a wide range of nuclear density. We 
then find a new, simple scaling relation for the changes of hadron masses in the medium: 

5ml 6MX^ 2 6 Mi 1 



5M% 5M% 3' 5M% 



etc. (1) 



where SM* = Mi — M*, with the effective hadron mass, M* {i = N,uj, p, ■ ■ ■). 

An outhne of the paper is as follows. In Sec. |T|, the idea of the QMC model is 
first reviewed. Then, the model is extended to include the effect of meson structure. 
In Sec. m, parameters in the model are first determined to reproduce the properties of 
infinite nuclear matter, and the hadron masses in the medium are then discussed. A new 
scaling relationship among them is also derived. The static properties of several closed- 
shell nuclei are studied in subsection [III.3i where we also show the changes of the masses 



of the nucleon, the mesons (cr, u and p) and the hyperons (A, S and H) in ^°Ca and ^°^Pb. 
The last section gives our conclusions. 

II THE QUARK-MESON COUPLING MODEL 

II. 1 Effect of nucleon structure in finite nuclei 

Let us suppose that a free nucleon (at the origin) consists of three light (u and d) quarks 
under a (Lorentz scalar) confinement potential, Vc- Then, the Dirac equation for the 
quark field, ipq, is given by 

[t'y- d - mg - Vc{r)]^g{r) = 0, (2) 

where rUg is the bare quark mass. 

Next we consider how Eq.(ED is modified when the nucleon is bound in static, uniformly 



distributed (iso-symmetric) nuclear matter. In the QMC model it is assumed that each 
quark feels scalar, V^, and vector, V^, potentials, which are generated by the surrounding 
nucleons, as well as the confinement potential (see also Ref.[Q). Since the typical distance 
between two nucleons around normal nuclear density (po = 0.15 fm~^) is surely larger 
than the typical size of the nucleon (the radius Rn is about 0.8 fm), the interaction 
(except for the short-range part) between the nucleons should be colour singlet; e.g., 
a meson-exchange potential. Therefore, this assumption seems appropriate when the 
baryon density, p_B, is not high. If we use the mean- field approximation for the meson 
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fields, Eq.(0) may be rewritten as 

[27 ■ 5 - K - K") - K(r) - 7oV7]V^,(r-0 = 0. (3) 

The potentials generated by the medium are constants because the matter distributes 
uniformly. As the nucleon is static, the time-derivative operator in the Dirac equation 
can be replaced by the quark energy, —ieg. By analogy with the procedure applied to 
the nucleon in QHD |]14|, if we introduce the effective quark mass by m* = rriq — V^, 



the Dirac equation, Eq.(§), can be rewritten in the same form as that in free space, with 
the mass m* and the energy — V^, instead of rriq and e^. In other words, the vector 
interaction has no effect on the nucleon structure except for an overall phase in the quark 
wave function, which gives a shift in the nucleon energy. This fact does not depend on 
how to choose the confinement potential, V^. Then, the nucleon energy (at rest), E^, in 



the medium is |jT9| 



E^ = M*^{V,^) + 3V,^, (4) 

where the effective nucleon mass, Mj^, depends on only the scalar potential in the medium. 

Now we extend this idea to finite nuclei. The solution of the general problem of a 
composite, quantum particle moving in background scalar and vector fields that vary 
with position is extremely difficult. One has, however, a chance to solve the particular 
problem of interest to us, namely light quarks confined in a nucleon which is itself bound 
in a finite nucleus, only because the nucleon motion is relatively slow and the quarks 
highly relativistic |jl6|. Thus the Born-Oppenheimer approximation, in which the nucleon 
internal structure has time to adjust to the local fields, is naturally suited to the problem. 
It is relatively easy to establish that the method should be reliable at the level of a few 
percent |]TB| . 



Even within the Born-Oppenheimer approximation, the nuclear surface gives rise to 
external fields that may vary appreciably across the finite size of the nucleon. Our ap- 
proach in Ref. was to start with a classical nucleon and to allow its internal structure 



to adjust to minimise the energy of three quarks in the ground-state of a system under 
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constant scalar and vector fields, with values equal to those at the centre of the nucleon. 



In Ref.|16], the MIT bag model was used to describe the nucleon structure. Blunden and 



Miller have also examined a relativistic oscillator model as an alternative model Of 
course, the major problem with the MIT bag (as with many other relativistic models of 
nucleon structure) is that it is difficult to boost. We therefore solve the bag equations in 
the instantaneous rest frame (IRF) of the nucleon - using a standard Lorentz transfor- 
mation to find the energy and momentum of the classical nucleon bag in the nuclear rest 
frame. Having solved the problem using the meson fields at the centre of the "nucleon" 
(which is a quasi-particle with nucleon quantum numbers), one can use perturbation the- 
ory to correct for the variation of the scalar and vector fields across the nucleon bag. In 
first order perturbation theory only the spatial components of the vector potential give 
a non- vanishing contribution. (Note that, although in the nuclear rest frame only the 
time component of the vector field is non-zero, in the IRF of the nucleon there are also 
non- vanishing spatial components.) This extra term is a correction to the spin-orbit force. 



As shown in Refs. ||T6|, 21 1, the basic result in the QMC model is that, in the scalar (a) 
and vector (uj) meson fields, the nucleon behaves essentially as a point-like particle with 
an effective mass M^, which depends on the position through only the a field, moving 
in a vector potential generated by the u meson, as mentioned near Eq.(^. Although we 



discussed the QMC model using the specific model, namely the bag model, in Ref.|16, 21 



the qualitative features we found are correct in any model in which the nucleon contains 
relativistic quarks and the (middle- and long-range) attractive and (short-range) repulsive 
N-N forces have Lorentz- scalar and vector characters, respectively. 

Let us suppose that the scalar and vector potentials in Eq.(^ are mediated by the a 
and u) mesons, and introduce their mean-field values, which now depend on position r, by 
V^{r) = glcr{r) and V^{f) = g^uj{r), respectively, where (g^) is the coupling constant 
of the quark-cr [u) meson. Furthermore, we shall add the isovector, vector meson, p, and 



the Coulomb field, A(r), to describe finite nuclei realistically |T^, Then, the effective 
Lagrangian density for finite nuclei, involving the quark degrees of freedom in the nucleon 
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and the (structureless) meson fields, in MFA would be given by 
Cqmc-i = ^[h ■ d - M^{a{r)) - g^uj{r)'yo 

- ^[(Va(f))2 + mla{r)'] + ^[(Vc.(r-))2 + mXf)^] 

+ ^[(V6(f)f + mJ6(fn + i(VA(f))^ (5) 

where ilj{f) and b{f) are respectively the nucleon and the p (the time component in the 
third direction of isospin) fields, m^, m^^ and irtp are respectively the (constant) masses 
of the a, uj and p mesons, g^) and gp are respectively the w-N and p-N coupling constants, 
which are related to the corresponding quark-tu, (7^, and quark-p, (7^, coupling constants 
as g^^ = 3g^ and gp = g^ [|16|, |2T[. We call this model the QMC-I model. If we define the 
field- dependent cr-N coupling constant, ga{<^), by 

M;^(a(r)) ^ - ^,(a(f))a(f), (6) 

where is the free nucleon mass, it is easy to compare with QHD [0]. gcr{o') will be 
discussed further below. 

The difference between QMC-I and QHD lies only in the coupling constant g^j, which 
depends on the scalar field in QMC-I while it is constant in QHD. (The relationship 
between QMC and QHD has been already clarified in Ref.|]19[. See also Ref.|2^.) However, 



this difference leads to a lot of favorable results, notably the nuclear compressibility. 



Tq, 18, IS, 21 1 . Detailed calculated properties of both infinite nuclear matter and finite 



nuclei can be found in Refs. [flBi pi| . 

Here we consider the nucleon mass in matter further. The nucleon mass is a function 
of the scalar field. Because the scalar field is small at low density the nucleon mass can 
be expanded in terms of cr as 

In the QMC model the interaction Hamiltonian between the nucleon and the a field at 
the quark level is given by Hint = —^qI I drip aipq, and the derivative of M^j with respect 



to a is 

m 

Here we have defined the quark-scalar density in the nucleon, S^la), which is itself a 
function of the scalar field, by Eq.(||). Because of a negative value of (-^), the nucleon 
mass decreases in matter at low density. 

Furthermore, we define the scalar-density ratio, Sj\f{a)/ 3^(0), to be Ciy{a) and the 
(T-N coupling constant at cr = to be g„ (i.e., g„ = g^i^a = 0)): 

Civ(a) = SN{(r)/SN{0) and g, = 3(7^^^(0). (9) 

Comparing with Eq.(^, we find that 

i-^j = -g.C^{<y) = -— [g^{a)a] , (10) 
and that the nucleon mass is 

M*^ = Mn - ga(T - ^g.C'^iO)a^ + ■ ■ ■ . (11) 

In general. Cat is a decreasing function because the quark in matter is more relativistic 
than in free space. Thus, C^(0) takes a negative value. If the nucleon were structureless 
Cm would not depend on the scalar field, that is, would be constant {C^ = 1). 
Therefore, only the first two terms in the right hand side of Eq.(^) remain, which is 
exactly the same as the equation for the effective nucleon mass in QHD. By taking the 
heavy- quark- mass limit in QMC we can reproduce the QHD results [ p!9|] . 

If the MIT bag model is adopted as the nucleon model, Sn is explicitly given by [|l^, |4 



^^^^^ " ?F(?F^TyTK^' ^ ^ 

where fi* = ^Jx*j^ + {R*j^m*y is the kinetic energy of the quark in units of 1/ and 
is the eigenvalue of the quark in the nucleon in matter. We denote the bag radius of the 
nucleon in free space (matter) by Rn (Rn)- actual numerical calculations we found 
that the scalar-density ratio. Cat (a), decreases linearly (to a very good approximation) 



with g^a |T6| , pT| . Then, it is very useful to have a simple parametrization for Cn'- 

CN{(y) = 1 - ttAT X {g„a), (13) 

with g„a in MeV (recall g^ = g^{a = 0)) and a^v ~ 9 x 10"^ (MeV^^) for = 5 MeV 
and i?Ar = 0.8 fm. This is quite accurate up to ~ 3po- 

As a practical matter, it is easy to solve Eq.(|^) for ga{cr) in the case where C{a) is 
linear in gcrCr, as in Eq.(|l^). Then one finds 



M* = - g^ 



1 / \ 

1 - ^(^af^J 



(14) 



so that the effective a-N coupling constant, ga{cr), decreases at half the rate of Cat (a). 
II. 2 Effect of meson structure 

In the previous section we have considered the effect of nucleon structure. It is however 
true that the mesons are also built of quarks and anti-quarks, and that they may change 
their properties in matter. 

To incorporate the effect of meson structure in the QMC model, we suppose that the 
vector mesons are again described by a relativistic quark model with common scalar and 
vector mean- fields |2^, like the nucleon (see Eq.(^). Then, again the effective vector- 



meson mass in matter, m*(f = uj,p), depends on only the scalar mean-field. 

However, for the scalar (cr) meson it may not be easy to describe it by a simple quark 
model (like a bag) because it couples strongly to the pseudoscalar (27r) channel, which 
requires a direct treatment of chiral symmetry in medium Since, according to the 



Nambu-Jona-Lasinio model ^ or the Walecka model [Q, one might expect the a- 
meson mass in medium, m*, to be less than the free one, we shall here parametrize it 
using a quadratic function of the scalar field: 

— =l-a,igaa) + b^{g^af, (15) 

with g^a in MeV, and we introduce two parameters, (in MeV""*^) and b^j (in MeV^^). 
(We will determine these parameters in the next section.) 
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Using these effective meson masses, we can find a new Lagrangian density for finite 
nuclei, which involves the structure effects of not only the nucleons but also the mesons, 
in the MFA: 

- r.^ e 

Cqmc-u = ^P[tj-d-M;,-g^uj{f)-fo-gp^b{r)^o--{l + T^)A{r)^o]^ 

- ^[(Va(f))2 + mfairr] + ^[(V^(f))2 + mlMm 

+ limf))' + mfbify] + ^(VA(f))^ (16) 

where the masses of the mesons and the nucleon depend on the scalar mean-fields. We 
call this model QMC-II. 

At low density the vector-meson mass can be again expanded in the same way as in 
the nucleon case (Eq.(|^): 



2 „ 1 



= - -g^ry/N(J - -gaTv/NC'^{0)o- , (17) 



where Sv{a) is the quark-scalar density in the vector meson, 

and Ct,(cr) = Sy{(r)/ 3^(0). In Eqs. (p^ and ([T8|), we introduce a correction factor, T^/j^, 
which is given by 3^(0)/ Sn{0), because the coupling constant, g^, is defined specifically 
for the nucleon by Eq.(|]). 



Ill NUMERICAL RESULTS 

In this section we will show our numerical results using the Lagrangian density of the 
QMC-II model that is, including self-consistently the density dependence of the meson 
masses. We have studied the QMC-I model, and have already shown the calculated 
properties of finite nuclei in Refs.[]TB|, pT| . 
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III.l Infinite nuclear matter 

For infinite nuclear matter we take tiie Fermi momenta for protons and neutrons to be 
{i = p or n). This is defined by pi = fc|n./(37r^), where pi is the density of protons or 
neutrons, and the total baryon density, ps, is then given by Pp + Pn- Let the constant 
mean-field values for the a, u and p fields be a, uj and 6, respectively. 

/^From the Lagrangian density Eq.(|l^), the total energy per nucleon, Etot/A, can be 
written (without the Coulomb force) 
2 ~ rkp. 



Etot/A 



Pb{2ti) 



E 

i=p,n ^ 



dk^Mf + P + 



+ -pB + 



9l 



2p 



B 



2m*2' 



SmfpB 



(19) 



where the value of the u) field is now determined by baryon number conservation as 
^ = QloPb , and the p-field value by the difference in proton and neutron densities, 



P3 = pp- pn, as 6 = gpp3/ i2mf) |20[. 

On the other hand, the scalar mean-field is given by a self-consistency condition (SCC): 



cr 



+ 



(27r)'^m*' 



*2 



E 

i=p,n 

(9m* ^ 
da I 



dk 



Mt 



+ 



Mf + fc2 
^ dm* 



'dM* 

~dw 



9ppl 



4:m*^m*^ \ da 



a 



dm* 



(20) 



Using Eqs.([TO|), (^) and (pD, Eq.(|OD can be rewritten 



a 



{2TT)^mf 



E C^(^) f 



dk 



M* 



i=p,n 



Mf + A:2 



2K{g,a)]a^ 



m 



*2 



'glpl 



^u/NCJa) + 



m 



Am*p^ 



^p/NCp{a) 



(21) 



Now we need a model for the structure of the hadrons involved. We use the MIT bag 



model in static, spherical cavity approximation [^. As in Ref. ||21||, the bag constant B 
and the parameter (which accounts for the sum of the cm. and gluon fluctuation 



corrections |jT6|) in the familiar form of the MIT bag model Lagrangian are fixed to 
reproduce the free nucleon mass (M^r = 939 MeV) under the condition that the hadron 
mass be stationary under variation of the free bag radius (i?7v in the case of the nucleon). 
Furthermore, to fit the free vector-meson masses, m^ = 783 MeV and m„ = 770 MeV, we 
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introduce new ^-parameters for them, z^^ and Zp. In the following we choose i?Ar = 0.8 fm 
and the free quark mass = 5 MeV. Variations of the quark mass and only lead to 
numerically small changes in the calculated resuhs |2l[. We then find that 5^/^ = 170.0 



MeV, = 3.295, z^ = 1.907 and Zp = 1.857. Thus, Cn is given by Eq.(|T2D, and is 
given by a similar form, with the kinetic energy of quark and the bag radius for the vector 
meson. We find that the bag model gives F^^p/Tv = 0.9996. Therefore, we may discard 
those correction factors in practical calculations. 

Next we must choose the two parameters in the parametrization for the cr-meson 
mass in matter (see Eq.(^)). In this paper, we consider three parameter sets: (A) 
a„ = 3.0 X 10"^ (MeV-i) and = 100 x 10"^ (MeV^^), (B) a„ = 5.0 x 10"^ (MeV-^) 
and 6^ = 50 X 10"^ (MeV-^), (C) = 7.5 x 10"^ (MeV~i) and = 100 x 10^^ (MeV-^). 
The parameter sets A, B and C give about 2%, 7% and 10% decreases of the a mass at 
saturation density, respectively. We will revisit this issue in the next subsection. 

Now we are in a position to determine the coupling constants, and g"^ are fixed to fit 
the binding energy (—15.7 MeV) at the saturation density (po = 0.15 fm~^) for symmetric 
nuclear matter. Furthermore, the p-meson coupling constant is used to reproduce the bulk 
symmetry energy, 35 MeV. We take = 550 MeV. The coupling constants and some 
calculated properties for matter are listed in Table ^ The last three columns show the 
relative changes (from their values at zero density) of the nucleon-bag radius {6R*p^/Rn), 
the lowest eigenvalue {6x*j^/xn) and the root-mean-square radius (rms radius) of the 
nucleon calculated using the quark wave function {Sr*/rq) at saturation density. We note 
that the nuclear compressibility is higher than that in QMC-I {K ~ 200 - 300 MeV) plj. 
However, it is still much lower than in QHD [Q. As in QMC-I, the bag radius of the 
nucleon shrinks a little, while its rms radius swells a little. On the other hand, because of 
the scalar field, the eigenvalue is reduced more than 10% (at po) from that in free space. 

The strength of the scalar mean-field, g^o", in medium is shown in Fig.|I[ At small 
density it is well approximated by a linear function of the density: 

^ 200 (MeV) ( — ) . (22) 



Po 
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Table 1: Coupling constants and calculated properties for symmetric nuclear matter at 
normal nuclear density {rUq — 5 MeV, — 0.8 fm and = 550 MeV). The effective 
nucleon mass, Mj^, and the nuclear compressibihty, K, are quoted in MeV. The bottom 
row is for QHD. 



type 


9l/^^ 


9l/^^ 


9l/^^ 




K 


SRf^/ Rn 






A 


3.84 


2.70 


5.54 


801 


325 


-0.01 


-0.11 


0.02 


B 


3.94 


3.17 


5.27 


781 


382 


-0.01 


-0.13 


0.02 


C 


3.84 


3.31 


5.18 


775 


433 


-0.02 


-0.14 


0.02 


QHD 


7.29 


10.8 


2.93 


522 


540 
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III. 2 A new scaling phenomenon for hadron masses in matter 

First, we show the dependence of the cx-meson mass on the nuclear density in Fig.^. Using 




0.5 1 1.5 2 2.5 3 

Pb^Po 



Figure 2: Effective cr-meson mass in symmetric nuclear matter. The curves are labelled 
as in Fig.l. 



Eqs.([T3|) and (P^, we find the a mass at low density is 



"'^-l-a, ^^^^ . (23) 



where = (0.06, 0.1, 0.15) for parameter set (A, B, C), respectively. 

The effective nucleon mass is shown in Fig.^. It decreases as the density goes up, and 
behaves like a constant at large density. At small density it is approximately given by 
using Eqs.(|l§ and (g^): 

^-^-l-0.2lH- (24) 



In Fig.^ the effective cj-meson mass is shown as a function of the density. (Since the 
difference between the effective u- and p-meson masses at the same density is very small, 
we show only one curve for both mesons in the figure.) As the density increases the vector- 
meson mass decreases (as several authors have previously noticed ||, ^ H, ^, |10|, |Tl], |12[) 
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Figure 3: Effective nucleoli mass in symmetric nuclear matter. The curves are labelled as 
in Fig.l. 



and seems to become ffat like the effective nucleon mass. Again, using Eqs.(|T^ and (p2D, 
the mass reduction can be well approximated by a linear form at small density: 

f^Ul-O.irH. (25) 

The reduction factor, 0.17, is consistent with other models which have been applied to 



the same problem |12 



In Fig. I we show the ratios of the quark-scalar density in medium to that in free space 
for the nucleon (Cat) and the u meson (C^). As pointed out previously, we can easily see 
that the ratio for the nucleon is well approximated by a linear function of (^o-cr. It is also 
true that the ratio for the vector meson can be well described by a similar, linear function 
of g^a: 

C^{a) = 1 - X (g^a). (26) 

We will see this parametrization again later. 

In the present model it is possible to calculate masses of other hadrons. In particular, 
there is considerable interest in studying the masses of hyperons in medium - eg. A, E 
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Figure 4: Effective (p- or) cj-meson mass in symmetric nuclear matter. The curves are 
labelled as in Fig.l. 

and S. For the hyperons themselves we again use the MIT bag model. We assume that 
the strange quark in the hyperon does not directly couple to the scalar field in MFA, as 
one would expect if the cr-meson represented a two-pion-exchange potential. It is also 
assumed that the addition of a single hyperon to nuclear matter of density ps does not 
alter the values of the scalar and vector mean-fields, namely, we take the local-density 
approximation to the hyperons The mass of the strange quark, m^, is taken to 

be mg = 250 MeV, and new ^-parameters in the mass formula are again introduced to 
reproduce the free hyperon masses: z/<^ = 3.131, z^, = 2.810, and z= = 2.860. Using those 
parameters, we have calculated the masses of A, E and S in symmetric nuclear matter. 
They are presented in Fig.^. As for the nucleon and the vector mesons, the effective mass 
of the hyperon is determined by only the scalar field. 

In general, we thus find that the effective hadron mass in medium is given by 




(27) 
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Figure 5: The ratios of the quark-scalar density in medium to that in free space for the 
nucleon (sohd curve) and the u) meson (dotted curve) - using parameter set B. 



where j stands for A^, u;, p, A, S, S, etc., no is the number of non-strange quarks in 
the hadron j, Tj/jv = S'j(0)/S'7v(0) with the quark-scalar density, Sj, in j, and the scalar 
density ratio, Cj(cr) = Sj{a)/Sj(fi). 

Using Eqs.(B^) and (pT]), we find that the hyperon masses at low density are given by 



Ml 
Ma 

My 



1 - 0.12 



and 



^1-0.11^ 
~ 1-0.05 I ^ 



(28) 
(29) 

.... , (30) 

VMh/ \PoJ 

where we take FA ^.s/Af = 1, because we find that the F factor for the hyperon is again 
quite close to unity (e.g. Fa/at = 1.0001, in our actual calculations). 

As seen in Fig.^ the linear approximation to the scalar-density ratio, Cj, is very 
convenient. We find that it is numerically relevant to not only the nucleon and the 
vector mesons but also the hyperons: 



Cj(cr) = l-aj X {g„a) 



(31) 
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Figure 6: The ratio of the hyperon mass in medium to that in free space. The dotted, 
sohd and dashed curves are respectively for the A, S, and S hyperons, using parameter 
set B. 



where aj is the slope parameter for the hadron j. We list them in Table |^. We should 



Table 2: Slope parameters for the hadrons (xlO ^ MeV ^). 



type 






ap 




as 




A 


9.01 


8.63 


8.59 


9.27 


9.52 


9.41 


B 


8.98 


8.63 


8.58 


9.29 


9.53 


9.43 


C 


8.97 


8.63 


8.58 


9.29 


9.53 


9.43 



note that the dependence of aj on the hadrons is quite weak, and it ranges around 8.6 ~ 
9.5 X 10"^ (MeV-^). 

If we ignore the weak dependence of aj on the hadrons and take ^j/n = 1 in Eg. (pTf) , 
the effective hadron mass can be rewritten in a quite simple form: 

no, 



1 - 2^9a(^^ 



(32) 
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where a ^ 9.0 x 10^^ (MeV~^). This mass formula can reproduce the hadron masses in 
matter quite well over a wide range of p^, up to ~ 3po- 

Since the scalar field is common to all hadrons, Eg . (|32D leads to a new, simple scaling 
relationship among the hadron masses: 

5ml\ ( 5MX\ f6M*\ 2 , f 6Mi\ 1 , , 

and — ^ - -, (33) 



\6M^J \6M^J \SM^J 3 3' 

where SM* = Mj — M*. The factors, | and |, in Eg .(^31) come from the ratio of the 
number of non-strange guarks in j to that in the nucleon. This means that the hadron 
mass is practically determined by only the number of non-strange guarks, which feel the 
common scalar field generated by surrounding nucleons in medium, and the strength of 



the scalar field . On the other hand, the change in the confinement mechanism due to 
the environment gives a small contribution to the above ratio. It would be very interesting 
to see whether this scaling relationship is correct in forthcoming experiments. 
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III. 3 Finite nuclei 

In this subsection we will show our results for some finite, closed shell nuclei. The La- 
grangian density, Eq.(|l^), leads to the following equations for finite nuclei: 

d? 2d 

—a{r) + -—a{r)-mfa{r) = -g„CNps{r) - m„mlg„[a^ - 2h„g„a{r)]a{rf 
dr^ r dr 

+ '^g,[mlV^,NCMrf + mlV^,^Cph{r)\ (34) 

(f 2d 

—uj{r) + -—uj{r)-m*Juj{r) = -gu^Psir), (35) 

^2bir) + l^^Kr)-mfb{r) = -|p3(r), (36) 

f-,A{r) + -^A{r) = -ep,{r), (37) 
dr^ r dr 



where 



Ps{r) = E^«W(|G.(r)p-|F.(r)p), (38) 

a 

occ 

PB{r) = j:Mr){\Ga{r)\' + \F^{r)f), (39) 



a 

occ 



P3{r) = T.d4r){-y-'^'{\G^{r)\' + \F4r)\'), (40) 

a 

occ -| 

Ppir) = j:Mr){ta + ^){\G^{r)\'+\F^{r)\'), (41) 

with djr) = {2ja + l)/47rr^, and 

d K 1 

—Gair) + -Gair) - [e„ - g^uj{r) - t^gpb{r) - {t^ + -)eA{r) + Mn 
dr r 2 

- g,{a{r))a{r)]F^{r) = Q, (42) 

-^F„(r) - -F,(r) + [e„ - g^uj{r) - t^g,b{r) - (t, + ^)eA(r) - 
dr r 2 

+ g^{a{r))a{r)]G^{r) = 0. (43) 

Here Ga{r)/r and Fa{r)/r are respectively the radial part of the upper and lower com- 
ponents of the solution to the Dirac equation for the nucleon: 

m = ( v., (44) 



■|F„(r)/r]*_ 



Km I 
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where C,ta is a two-component isospinor and is a spin spherical harmonic |^9|] {a la- 
belhng the quantum numbers and being the energy). Then, the normahzation condition 
is 

|dr(|G,(r)p + |F„(r)p) = l. (45) 

As usual, K specifies the angular quantum numbers and the eigenvalue of the isospin 
operator Practically, m*, m* and Cj are respectively given by Eqs.(|l^), (p7|) and 



(pip, and gfj{a{r)) is 



1 - f{9Mr)) 



(46) 



The total energy of the system is then given by 

OCC p 

Etot = JZC^Ja + l)e„ -^J dr [-g,D{a{r))a{ 



+ gu^^{.r)pB{r) + ]^gph{r)p^{r) + eA{r)pp{r)], (47) 



where 



D{a{r)) = CNps{r) + m„ml[a^ -2b„g„a{r)\a{rf 
2 
3 



There are seven parameters to be determined: g^, g^i Qp, e, rria-, and nip. As in 
the case of infinite matter we take the experimental values: m^^ = 783 MeV, rup = 770 
MeV and e^/47r = 1/137.036. The coupling constants g^, gui and gp are fixed to describe 
the nuclear matter properties and the bulk symmetry energy per baryon of 35 MeV (see 
Table I). 

The (T-meson mass however determines the range of the attractive interaction and 
changes in mg- affect the nuclear-surface slope and its thickness. Therefore, as in the 
paper of Horowitz and Serot |^, we adjust to fit the measured rms charge radius of 
^°Ca, rc/i('^°Ca) = 3.48 fm (Notice that variations of mg. at fixed {ga/^a) have no 



effect on the infinite nuclear matter properties ||2T[].) We summarize the parameters in 
Table |3|. 



21 



Table 3: Model parameters for finite nuclei (for rriq = 5 MeV and Rn = 0.8 fm). 



type 








m^(MeV) 


A 


1.67 


2.70 


5.54 


363 


B 


2.01 


3.17 


5.27 


393 


C 


2.19 


3.31 


5.18 


416 




Figure 7: Scalar and vector strength for ^°Ca and ^°^Pb (for type B). 



Equations (|35D to (|^) give a set of coupled non-linear differential equations, which 
may be solved by a standard iteration procedure ||3ll . 



In Fig.^we first show the calculated strength of the a and uj fields in ^°Ca and ^°^Pb. 
Next we show calculated charge density distributions, pch, of ^°Ca and ^^'^Pb in comparison 
with those of the experimental data in Figs.^ and ^. To see the difference among the 
results from the three parametrizations of m* (A, B and C), in Fig.^we present only the 
interior part of Pc/i (^°Ca). As in Ref. ||21||, we have used a convolution of the point-proton 
density, which is given by solving Eqs.(|35|) ~ (^51), with the proton charge distribution 
to calculate pch- For ^°Ca the QMC-11 model with parameter sets A and B give similar 
charge distributions to those in QMC-I, while the result of QMC-II with parameter set C 
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Figure 8: Charge density distribution for "^^Ca compared with the experimental data. 

is cfoser to that in QHD. From Fig.^ we see that the present model also yields a charge 
distribution for ^°'^Pb which is similar to those calculated using QMC-I or QHD. 

In Table ^, the calculated spectrum of "^^Ca is presented. Because of the relatively 
smaller scalar and vector fields in the present model than in QHD, the spin-orbit splittings 
are smaller (in this respect the model is very similar to QMC-1). We should note that there 



is a strong correlation between the effective nucleon mass and the spin-orbit force |21 



The problem concerning the spin-orbit force in the QMC model has been studied in 



Refs. jig, |21| , P^ , |23|| . It remains to be seen whether the higher order corrections, as 



studied by Wallace et al. will help to resolve it. 

Table ^ gives a summary of the calculated binding energy per nucleon {E/A)^ rms 
charge radii and the difference between nuclear rms radii for neutrons and protons (r„— r^), 
for several closed-shell nuclei. While there are still some discrepancies between the results 
and data, the present model provides reasonable results. In particular, as in QMC-I, it 
reproduces the rms charge radii for medium and heavy nuclei quite well. 

In Figs.|^ and [Tl| we present the changes of the nucleon, a and uj meson masses in 
^°Ca and ^°^Pb, respectively. The interior density of ^''Ca is much higher than p^, while 
that in ^''^Pb is quite close to po- Accordingly, in the interior the effective hadron masses 
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Table 4: Calculated proton and neutron spectra of ^'^Ca (for type B) compared with 
QMC-I and the experimental data {rriq = 5 MeV and Rn = 0.8 fm). Here, I and II 
denote, respectively, QMC-I and QMC-II. All energies are in MeV. 







neutron 




proton 


Shell 


I 


II 


Expt. 


I 


II 


Expt. 




43.1 


41.1 


51.9 


35.2 


33.2 


50±10 


lP3/2 


31.4 


30.0 


36.6 


23.8 


22.3 


34±6 


lPl/2 


30.2 


29.0 


34.5 


22.5 


21.4 


34±6 


14/2 


19.1 


18.0 


21.6 


11.7 


10.6 


15.5 


2Si/2 


15.8 


14.7 


18.9 


8.5 


7.4 


10.9 


14/2 


17.0 


16.4 


18.4 


9.7 


9.0 


8.3 



Table 5: Binding energy per nucleon, —E/A (in MeV), rms charge radius Vch (in fm) and 
the difference between r„ and (in fm) for type B, iriq = 5 MeV and Rb = 0.8 fm. I 
and II denote, respectively QMC-I and QMC-II. (* fit) 





-E/A 


rch 




model 


I 


II 


Expt. 


I 


II 


Expt. 


I 


II 


Expt. 


160 


5.84 


5.11 


7.98 


2.79 


2.77 


2.73 


-0.03 


-0.03 


0.0 


^oCa 


7.36 


6.54 


8.45 


3.48* 


3.48* 


3.48 


-0.05 


-0.05 


0.05±0.05 


^«Ca 


7.26 


6.27 


8.57 


3.52 


3.53 


3.47 


0.23 


0.24 


0.2±0.05 


902r 


7.79 


6.99 


8.66 


4.27 


4.28 


4.27 


0.11 


0.12 


0.05±0.1 


208 


7.25 


6.52 


7.86 


5.49 


5.49 


5.50 


0.26 


0.27 


0.16±0.05 



24 



0.08 




Figure 9: Same as Fig.8 but for ^''^Pb. The parameter set B is used. 

in ^°Ca are smaller than in ^°^Pb. We can also see this in Fig.|^, where the strength of the 
scalar field in the interior part of '^°Ca is stronger than in ^''^Pb. 

Using the local-density approximation and Eg. (|32D , it is possible to calculate the 
changes of the hyperon (A, S and H) masses in '^''Ca and ^°^Pb, which are respectively 
illustrated in Figs.^ and 0. Our quantitative calculations for the changes of the hyperon 
masses in finite nuclei may be quite important in forthcoming experiments concerning 



hypernuclei [2S 
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r(fm) 

Figure 10: Changes of the nucleon, a and oj meson masses in ^°Ca. The nuclear baryon 
density is also illustrated (solid curve). The right (left) scale is for the effective mass (the 
baryon density). The parameter set B is used. 
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Figure 12: Changes of the hyperon (A, E and S) masses in ^'^Ca. The sohd curve is for 
the nuclear baryon density. The right (left) scale is for the effective mass (the baryon 
density). The parameter set B is used. 




IV CONCLUSION 



We have extended the quark-meson couphng (QMC) model to include quark degrees 
of freedom within the scalar and vector mesons, as well as in the nucleons, and have 
investigated the density dependence of hadron masses in nuclear medium. As several 
authors have suggested [|, 0, ||, ^ 0' 0' HI' 1131 ' hadron mass is reduced because of 



the scalar mean-field in medium. Our results are quite consistent with the other models. 
In the present model the hadron mass can be related to the number of non-strange quarks 
and the strength of the scalar mean-field (see Eg . (|32|) ) . We have found a new, simple 
formula to describe the hadron masses in the medium, and this led to a new scaling 
relationship among them (see Eq.(p3D). Furthermore, we have calculated the changes of 
not only the nucleon, a, uj and p masses but also the hyperon (A, S and S) masses in 
finite nuclei. We should note that the origins of the mass reduction in QMC and QHD 
are completely different [|1^. It would be very interesting to compare our results with 
forthcoming experiments on hypernuclei. 

By applying this extended QMC model to finite nuclei, we have studied the properties 
of some static, closed shell nuclei. Our (self-consistent) calculations reproduce well the 
observed static properties of nuclei such as the charge density distributions. In the present 
model, there are, however, still some discrepancies in energy spectra of nuclei, in partic- 
ular, the spin-orbit splittings. To overcome this defect, we have discussed one possible 
way, in which a constituent quark mass (~ 300 MeV) is adopted, in Refs. |I6|, |2T|. As an 



alternative, Jin and Jennings p3[ and Blunden and Miller p2[ have proposed variations 
of the bag constant and z parameter in medium, which have been suggested by the fact 
that quarks are partially deconfined in matter. To help settle this problem, one should 
perhaps consider the change of the vacuum properties in the medium p5| . 

Our Lagrangian density, Eq. ([l6|) , provides a lot of effective coupling terms among 
the meson fields because the mesons have structure (cf. Ref. ||33||). In particular, the 
Lagrangian automatically offers self-coupling terms (or non-linear terms) with respect to 
the a field. Using Eq.(^), the Lagrangian density gives the non- linear a terms (up to 

28 



0{a^)) as: 



- -\rnla^ + g„a^m'y ~]^gl{al + 2K)mla\ (49) 

On the other hand, in nuclear physics, QHD with non-hnear a terms has been exten- 
sively used in MFA to describe realistic nuclei [0] . The most popular parametrizations are 



called NLl, NL2 |35| and NL-SH ^6[, and the non-linear terms in those parametrizations 
are given as 



i^QHD = -l^W + 192^' + ^93^', (50) 



where g2 and g^ take respectively a positive (negative) [positive] and positive (negative) 
[positive] values in NLl (NL2) [NL-SH]. Since the non-linear a terms provide the self- 
energy of a meson, it changes the a mass in matter. Comparing Eq.(pOD with Eg . (^9]) , 
we can see that the effective a mass in NL2 increases at low nuclear density while the a 
mass decreases in NLl and NL-SH in MFA. 

However, from the point of view of a field theory, like the Nambu-Jona-Lasinio model, 
an increase of the a mass in the medium seems unlikely (We should note that the 

values of g2 in those parametrizations are small compared with the corresponding one in 
Eq.(^9]).) Furthermore, from the point of view of field theory, g^ in Eg. (pO]) should be neg- 
ative because the vacuum must be stable |3^. Therefore, we can conclude that one would 
expect to find (72 > and (73 < in Eg. (|50|) . Unfortunately, the above three parametriza- 
tions used in nuclear physics do not satisfy the condition, while our Lagrangian, Eg. (^9]) , 
does. It will be very interesting to explore the connection between various coupling 
strengths found empirically in earlier work and those found in our approach. 

Finally, we would like to give some caveats concerning the present calculation. The 
basic idea of the model is that the mesons are locally coupled to the guarks. Therefore, in 
the present model the effect of short-range correlations among the guarks, which would 
be associated with overlap of the hadrons, are completely neglected. At very high density 
these would be expected to dominate and the present model must eventually break down 
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there (probably beyond ~ Sps/po)- Furthermore, the pionic cloud of the hadron 
should be considered explicitly in any truly quantitative study of hadron properties in 
medium. We note that subtleties such as scalar- vector mixing in medium and the splitting 
between longitudinal and transverse masses of the vector mesons [0 have been ignored 
in the present mean-field study. Although the former appears to be quite small in QHD 
the latter will certainly be important in any attempt to actually measure the mass shift. 

This work was supported by the Australian Research Council. 
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